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ABSTRACT 

The color-magnitude diagram (CMD) of the intermediate-age Large Magellanic Cloud 
(LMC) star cluster NGC 2154 and its adjacent held, has been analyzed using Padova 
stellar models to determine the cluster's fundamental parameters and its Star For- 
mation History (SFH). Deep BR CCD photometry, together with synthetic CMDs 
and Integrated Luminosity Functions (ILFs), has allowed us to infer that the cluster 
experienced an extended star formation period of about 1.2 Gyrs, which began ap- 
proximately 2.3 Gyr and ended 1.1 Gyr ago. The physical reality of such a prolonged 
period of star formation is however questionable, and could be the result of inadequa- 
cies in the stellar evolutionary tracks themselves. A substantial fraction of binaries 
(70%) seems to exist in NGC 2154. 

Key words: color-magnitude diagrams - galaxies individual (Large Magellanic 
Cloud) -galaxies:star clusters - star clusters: individual (NGC2154)- stars: evolution 



X 



1 INTRODUCTION 

The star cluster system of the Magellanic Clouds (MCs) dif- 
fers significantly from that of the Milky Way (and also from 
one another), differences which are commonly attributed to 
a different chemical and dynamical evolution. Furthermore, 
MC clusters exhibit a broad range of properties in contrast 
to our galaxy, thus representing a more ample range of stel- 
lar populations than those represented by Galactic clusters. 
For the above reasons, MCs clusters have become a chal- 
lenging domain for stellar and galactic evolutionary models, 
and are routinely used as an observational workbench to 
address these issues (see e.g. Barmina et al. 2002, Bertelli 
et al. 2003, Woo et al. 2003). A very specific case is that 
of the intermediate-age, metal-poor populations, which are 
conspicuous in the MC, yet rather poorly represented in our 
galaxy. 

One of many examples of intermediate-age clusters in 
the LMC is NGC 2154. Although this cluster is morpho- 
logically globular, it is considered to be of intermediate age 



(SMB type V; Searle et al. 1980). Persson et al. (1983) were 
able to constrain its age to the 1-3 Gyr range, but a detailed 
study of its basic parameters (requiring high-precision deep 
photometry) was still lacking. 

Here we present deep CCD photometry, reaching 
B,R~25, of NGC 2154 and its adjacent LMC field, which 
has allowed for an unprecedented study of the cluster, and 
a first-ever study of the SFH of this LMC field. The present 
work is one result of a more comprehensive study of the 
MCs, which includes the study of their SFH (Noel et al. 
2006; submitted to AJ), and the determination of their ab- 
solute proper motions with respect to background QSOs (see 
e.g. Pedreros et al. 2006). One of the LMC QSO fields se- 
lected for the proper-motion work by chance included the 
neglected LMC cluster NGC 2154, which gave us the pos- 
sibility to observe this cluster -and its surrounding field- 
routinely during our four-year (2001-4) campaign, with the 
(additional) motivation of determining not only its funda- 
mental parameters, but also the SFH of the field (which will 
be the subject of a forthcoming paper; Girardi et al. 2007). 



* email: gbaume@fcaglp.fcaglp.unlp.edu.ar 
(GB), giovanni.carraro@unipd.it (GC), costa@das.uchile.cl (EC), 
rmendez@das.uchile cl (RAM), leo.girardi@oapd.inaf.it (LG) 



The layout of the paper is as follows. Sect. 2 and 3 de- 
scribe the observations and the reduction strategy. In Sect. 4 
we study the cluster structure and derive an estimate of its 
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Figure 1. R band image of NGC 2154 and the adjacent LMC stellar field. North is up, East to the left. The size of the field is 8f85 X 8.' 85 



radius. Sect. 5 and 6 deal with the CMDs and describe the 
derivation of the cluster fundamental parameters. In Sect. 7 
we summarize our conclusions. 



2 OBSERVATIONS 

B(R)kc images of the NGC 2154 region in the LMC were 
acquired with a 24/i pixels Tektronix 2048x2048 detector 
attached to the Cassegrain focus of the du Pont 2.5-meter 
telescope (C100) at Las Campanas Observatory, Chile. Gain 
and read noise were 3 e-/ADU and 7 e-, respectively. This 



set-up provides direct imaging over a field of 8f85 x 8.'85 
with a scale of O'.'259/pix. This relatively large field of view 
allowed us to include the cluster and a good sampling of the 
LMC field in all frames. The field covered by the observa- 
tions is shown in Fig. 1. 

Details on the available frames and their corresponding 
exposure times are listed in Table 1. The B and R band- 
passes were selected in order to satisfy both the needs of 
the SFH and astrometry programs (for this latter it was 
mandatory to obtain R-band images). Typical seeing was 
about 0'.'9. 

All frames were pre-processed in a standard way using 
the IRAFt package CCDRED. For this purpose, zero expo- 
sures and sky flats were taken every night. 

3 THE PHOTOMETRY 

3.1 Standard star photometry 

Our instrumental photometric system was defined by the 
use of the Harris UBVRI filter set, which constitutes the 
default option at the C100 for broad-band photometry on 
the standard Johnson-Kron- Cousins system. On photomet- 
ric nights, standard star areas from the catalog of Landolt 
(1992) were observed multiple times to determine the trans- 
formation equations relating our instrumental (b, r) magni- 
tudes to the standard (B, Kkc) system. To determine at- 
mospheric extinction optimally, a few of them were followed 
each night up to about 2.0 air-masses. Besides, the stan- 
dard star fields were selected to provide a wide color 
coverage, being —0.4 < (B — R) < 3.6 (see Fig. 2). Aper- 
ture photometry was carried out for all the standard stars 
using the IRAF DAOPHOT/PHOTCAL package. 

To put our observations into the standard system, we 
used transformation equations of the form: 

6 = B + 61 + 62 * X + 63 * (B - R) 
r = R + rl + r2 * X + r3 * (B - R) 

In these equations 6, r are the aperture magnitude al- 
ready normalized to 1 sec, and X is the airmass. We did not 
include second-order color terms because they turned out 
to be negligible in comparison to their uncertainties. The 
values of the transformation coefficients are listed in Table 
2. The night-to-night variation of the coefficients turned out 
to be very small (~ 0.001), so we adopted the average val- 
ues over all nights. The residuals resulting from the 
fits are showed in Fig. 2, and the global rms of the 
calibration was 0.011 mag for both B and R niters. 

3.2 Cluster and LMC field photometry 

Here we follow the procedure outlined in Baume et al. 
(2004). We first averaged images taken the same night, and 

t IRAF is distributed by NOAO, which are operated by AURA 
under cooperative agreement with the NSF. 



The star cluster NGC 2154 

Table 1. Log-book of observations for NGC 2154. 



Date 


Airmass 


Filter 


Exp. Time 








[sec.xN] 


15-10-01 


1.28 


R 


60x1 




1.28 


R 


300x3 


is in 01 


1 97 


/ 1 


/inn v °. 


17-10-01 


1.28 


B 


400x1 




1.27 


R 


400x3 


18-10-01 


1.28 


R 


400x3 


08-10-02 


1.51 


R 


60x 1 




1.49 


R 


400x9 




1.30 


B 


400x1 




1.29 


B 


600x4 


09-10-02 


1.51 


R 


60x2 




1.31 


B 


60x1 




1.30 


B 


600x1 




1.28 


R 


300x3 




1.28 


R 


400x1 


10-10-02 


1.29 


B 


60x1 




1.29 


B 


600x1 




1.28 


R 


400x1 


11-10-02 


1.37 


B 


60x1 




1.32 


B 


600x4 




1.30 


R 


400x2 




1.28 


R 


300x4 


12-10-02 


1.37 


B 


60x1 




1.37 


B 


600x3 




1.28 


R 


60 x 1 




1.28 


R 


OUU A O 


20-10-03 


1.29 


R 


60x1 




1.29 


f> 


200x4 




1.27 




800 x 3 


21-10-03 


1.29 


R 


60x1 




1.29 


jj 


OUU A O 




1.28 


R 


300x4 


22-10-03 


1.29 


R 


600x1 




1.27 


R 


300x4 


24-10-03 


1.31 


R 


60x1 




1.31 


R 


600x2 




1.29 


R 


300x6 


04-11-04 


1.30 


R 


60x1 




1.29 


R 


600x6 




1.28 


B 


800x3 


05-11-04 


1.31 


R 


5x3 




1.30 


R 


30x3 




1.30 


R 


300x2 




1.29 


B 


10x4 




1.28 


B 


120x3 




1.28 


B 


600x2 




1.27 


R 


60x2 




1.27 


R 


300x5 




1.27 


B 


600x6 



N indicates the number of frames obtained. 
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Table 2. Coefficients of the calibration equations 



Table 3. Completeness analysis results in the B band. 



run 2001 
(four nights) 



bi 


+0.734 ±0.003 


hi 


+0.218 ±0.002 




-0.044 + 0.001 


n 


+0.434 + 0.003 




+0.083 ± 0.002 


rz 


-0.001 + 0.001 



with the same exposure time and filter, in order to get 
a higher signal-to-noise ratio for the faint stars, and also 
to clean the images from cosmic rays. Then, instrumental 
magnitudes and (X,Y) centroids of all stars in each frame 
were derived by means of profile-fitting photometry with 
the DAOPHOT package, using the Point Spread Function 
(PSF) method (Stetson 1987); and all instrumental mag- 
nitudes from different nights were combined and 
carried to same system (2001 reference) by using 
DAOM ASTER (Stetson 1992). Stellar magnitudes 
in the standard system were obtained then by using 
the transformations indicated in section 3.1. 
This resulted in a photometric catalog consisting of about 
20000 stars 

In Fig. 3 we present the photometric errors trends (from 
DAOPHOT and DAO MASTER) as a function of the B mag- 
nitude. Down to B = 22, both the B band and R band errors 
remain lower than 0.025 mag. From Fig. 2 it can be easily 
seen that the main source of error originates in the B-band 
observations. It is this filter that determines the deepness of 
our photometry, a point that is relevant to the completeness 
analysis (see Sect. 3.3). 

3.3 Photometric completeness 

For all comparisons described in the next sections, complete- 
ness of the observed star counts is a relevant issue. Complete- 
ness corrections were determined by means of artificial-star 
experiments on our data (see Carraro et al. 2005). Basi- 
cally, we created several artificial images by adding 
in random positions a total of 40000 artificial stars 
to our true images. They were distributed with a 
uniform probability distribution with the same color 
and luminosity of the real sample . In order to avoid 
the creation of overcrowding, in each experiment we added 
the equivalent to only 15of the original number of stars. 

Given that in general the B band images are shallower 
than those in the R band, we have adopted as completeness 
factor that estimated for B. This factor is defined as the 
ratio between the number of artificial stars recovered and 
the number of artificial stars added. Computed values of 
the completeness factor for different B magnitude bins are 
listed in Table 3, both for the cluster (r j 400 pix), and for 
a representative comparison field (see Sect. 6.1). It must 
be noticed that due to the inherent nature of a very 
compact star cluster (see Figs. 1 and 4), more than 
half of the stars would occupy less than half of the 
volume. As a consequence, the completeness frac- 
tions for the cluster stars are likely to be a bit over- 



AB 


NGC 2154 
(r<400pix) 


Field 


19.5-20.0 


100.0% 


100.0% 


20.0-20.5 


92.7% 


100.0% 


20.5-21.0 


75.1% 


100.0% 


21.0-21.5 


57.4% 


100.0% 


21.5-22.0 


56.8% 


100.0% 


22.0-22.5 


56.1% 


100.0% 


22.5-23.0 


55.0% 


100.0% 


23.0-23.5 


53.9% 


82.7% 


23.5-24.0 


41.7% 


61.1% 


24.0-24.5 


29.5% 


62.0% 


24.5-25.0 


32.4% 


65.6% 


25.0-25.5 


35.2% 


82.5% 


25.5-26.0 


36.3% 


49.6% 
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Figure 2. Trend of residuals of the standard stars calibration. 



estimated. On the other hand, the use of a lower 
radius for the cluster region would have the disad- 
vantage to imply larger uncertainties in the com- 
pleteness factors. 



4 STAR COUNTS 
4.1 Cluster radius 

In order to study the cluster structure, as a first step we es- 
timated the position of the cluster center by determining the 




Figure 3. Photometric errors in B and R given by DAOPHOT 
and DAOMASTER, as a function of the B magnitude. 
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highest peak in the stellar density, which was done by visual 
inspection of the our images. This peak was found at: X — 
515; Y = 1725, which corresponds to 012000 = 5 h 57 m 38 s .2; 
#2000 = — 67°15'40".7, coordinates which are similar to those 
given in the SIMBAD database. 

The next step was to compute the cluster's size, which 
was done constructing radial profiles by two methods: the 
radial stellar density profile and the the radial flux profile 
methods: 

a) In the first, stars are counted in a number of succes- 
sive rings, 30" wide, concentric around the adopted cluster 
center, and then divided by their respective areas. Because 
our data does not permit complete annuli beyond 515 pix, 
we have assumed that the measurable annuli portions are 
still representative of the field stellar populations around 
the cluster. The density profiles obtained, down to two dif- 
ferent B limit magnitude limits (20 and 21) are shown in 
Fig. 4a. 

b) In the second 
method, the flux (—2.5 log [ADU s / area]) within concentric 
annuli 10 pixels wide (2'.'59) is measured directly over the 
(B band) cluster image. The resulting profile is presented 
in Fig. 4b. A measure of the cluster's radius is obtained by 
fitting a model from Elson, Fall and Freeman (EFF; see El- 
son et al. 1987), appropriate for LMC clusters (Mackey & 
Gilmore 2003). The expression used was: 



H = Ho[l + (r/a) 2 ] 



- 7 /2 . 



where r is the distance from the adopted cluster center, fj, a 
the central surface brightness, a is a measure of the core 
radius, 7 the power law slope at large radii, and, finally, is 
the field surface brightness. The computed parameters are 
given in Fig. 4b. 



4.2 Observed Luminosity Function 

We have constructed B-band Luminosity Functions (LFs) of 
the cluster region (r < 400 pix) and of a similarly size field 
region (centered at X = 1500, Y = 500). The counts pre- 
sented were corrected using the completeness factors given 
in Table 3 (see Sect. 3.3). The results are plotted in Fig. 5a. 
In Figure 5b we present the pure cluster LF obtained by 
subtracting the field region LF from the cluster region LF. 

Finally, we separated the data (both in the cluster re- 
gion and in the field region) in two sets in order to isolate the 
red clump (RC) stars: those with B < 21.5 and B-R > 1.0, 
from the rest of the data (the main sequence region, MS). 
The resulting, field subtracted, LFs for the RC and MS re- 
gions are given in Fig. 5c and Fig. 5d, respectively. 



5 COLOR-MAGNITUDE DIAGRAMS 

In Fig. 6 we present the CMD of all stars measured in the 
complete field shown in Fig. 1. In Fig. 7 we present the 
CMD of the region centered at X=515, Y=1725 and hav- 
ing r < 400 pix (the cluster region, see Sect. 4). Although 
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Figure 4. a) Radial density profiles for NGC 2154. Numbers in- 
dicate the limit B magnitudes for each case, b) Radial flux profile 
for NGC 2154 (grey circles) and the EFF model fit (solid curve) 
together with the computed fitted parameters (c is an arbitrary 
constant). 



well centered in the cluster, this CMD is clearly contami- 
nated by LMC field stars (it is interesting to note that, as it 
is, this CMD closely resembles that of the intermediate-age 
LMC cluster NGC 2173 studied by Gallart et al. (2003)). 

To obtain a cleaner CMD, we applied the statistical de- 
contamination method described in Vallenari et al. (1992) 
and Gallart et al. (2003). In this procedure, a statistical 
subtraction of field stars is carried out making a star-by- 
star comparison between selected reference field regions and 
the cluster region. Briefly, for any given star in the reference 
field regions we look for the most similar (in color and mag- 
nitude) star in the cluster region, and remove it from the 
cluster's CMD. It should be noted that the procedure takes 
into account the different completeness level of the cluster 
and the field. 

For the above purpose, we selected three reference field 
regions having the same area of the cluster region, which 
we call (see Fig. 8) Field#l (centered at X= 1500,Y=500), 
Field#2 (centered at X=1500,Y=1500) and Field#3 (cen- 
tered at X=500,Y=500). They were chosen at proper dis- 
tances from the cluster center, in order to avoid the presence 
of cluster stars in them. The results are shown in Fig. 8. The 
three upper panels show the CMDs of the reference fields, 
whereas the middle and lower panels show the correspond- 
ing CMDs of the subtracted stars, and the corresponding 



6 Baume et al. 



16 
IB 

24 
IB 



j_ i | \ r i^^r"- fff4f 



I I I I I I I I I I ' I I I I 
Field J 





Li | I I 1 1 | I M I | 
■ .Field 




| i i i i-l i i i ' I i. 
Field 3 




1 2 
B-E 



Figure 8. Field star decontamination procedure. The upper panel shows the CMDs of three different field regions. The middle panel 
shows the corresponding CMDs of the subtracted stars. In the lower panel the corresponding clean cluster CMDs are shown. See text 
for details. 



clean cluster CMDs, respectively. 

Probably due to peculiarities in the distribution of field 
stars across the cluster region, a perfect decontamination 
was however not possible. Inspection of Fig. 8 shows that 
clean cluster regions present several groups of stars that are 
more numerous than in the reference field regions, despite 
the fact that all regions have the same area. In the clean 
CMDs some field stars still remain above the MS, to the left 
of the MS and everywhere in the red clump region. This ef- 
fect surely results from the statistically low number of stars 
in those regions of the CMD. 

Nonetheless, the procedure was very effective and 
helped us to improve the shape of the turnoff region (TO), 
and to remove several field red giant Branch (RGB) and 
horizontal branch stars. A careful inspection of Fig. 8 led us 
to adopt as the clean NGC 2154 CMD, that obtained us- 



ing Field #3 for the statistical decontamination (lower right 
panel). In this case, the subtracted field is the closest to the 
corresponding original field, and the number of field stars 
still in the cluster region is significantly lower than in the 
other two cases. 



6 COMPARISON WITH STELLAR MODELS 

In this section we derive estimates of NGC 2154's funda- 
mental parameters by comparing its CMD with theoretical 
stellar evolutionary models from the Padova library of stel- 
lar tracks and isochrones (Girardi et al. 2000). These mod- 
els have already been used in the past to study LMC star 
clusters with satisfactory results (e.g. Barmina et al. 2002, 
Bertelli et al. 2003). We first summarize previous work on 
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Figure 5. a) LFs corrected by completeness factors. The grey 
histogram corresponds to the 'cluster area', whereas the white 
histogram to the 'field area', b) LF corrected by field contamina- 
tion, c) and d) Idem as panel b) but only for RC and MS regions 
respectively, (see Sect. 4.2). 
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Figure 6. NGC 2154 CMD for all the measured stars. 
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Figure 7. NGC 2154 CMD for the stars lying within 400 pixel 
from the cluster center. The solid line encloses the area where we 
looked for evolved stars. 



NGC 2154 and what we know from the literature of its basic 
parameters. 



6.1 Metallicity 

Bica et al. (1986) used H/3 and G-band photometry to es- 
timate that the metallicity of NGC 2154 is Z = 0.006. 
NGC 2154 was later observed by Olszewski et al. (1991) 
as part of a spectroscopic survey of giant stars in LMC star 
clusters, who derived a metallicity measuring the pseudo- 
equivalent width of three Ca lines. Two stars were mea- 
sured with this purpose, giving an estimate of the cluster's 
metallicity of [Fe/H]=-0.56±0.20. Adopting the Carraro et 
al. (1999) relation, this value corresponds to Z = 0.005. We 
shall adopt this estimate of the metallicity throughout this 
work. 



6.2 Reddening and Distance modulus 

While the distance modulus to the LMC is known with rea- 
sonable precision (Westerlund 1997, (V — My) = 18.5), no 
estimates of the reddening in the direction of NGC 2154 
are available. To complicate matters, the reddening across 
the galaxy is known to be highly variable (Oestreicher & 
Schmidt-Kaler 1996; Zaritsky et al. 2004). As for the Galac- 
tic extinction law, here we shall use Rv—3.1 (Rieke & Lebof- 
sky 1985). 



6.3 Isochrone fitting 

In Fig. 9 we have superposed three Z=0.005 ([Fe/H]= -0.60) 
isochrones, taken from the Girardi et al. (2000) database, to 
our adopted clean NGC 2154 CMD. It should be noted that 
this metallicity is not directly available, so it has been inter- 
polated from different metallicity datasets. We have selected 
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Figure 9. CMD of NGC 2154 decontaminated from field stars 
(see see Sect. 5). The three Z = 0.005 isochrones superposed are 
from the Girardi et al. (2000) database (see Sect. 6.3). 

these isochrones because they provide a good fit to the MS 
and the TO region, and also to the magnitude and color of 
the RC. They have been shifted by the reddening and dis- 
tance modulus indicated in the figure. The values presented 
for these parameters imply a corrected distance modulus of 
(B — Mb) in the range 18.45 to 18.50, in good agreement 
with the widely adopted distance modulus given by Wester- 
lund (1997). 

Quick examination of Fig. 9 shows that: 

• The MS is significantly wide. This may be due to a com- 
bination of observational errors, binarity, differential redden- 
ing, and age spread; 

• The RGB clump is wide in color, and slightly tilted. 
This again can be ascribed to differential reddening and bi- 
narity; 

• The stars above the TO, and outside the MS edges are 
mostly interlopers belonging to the LMC field in front of the 
cluster; 

• On the other hand, the reddest field stars are probably 
members of the Milky Way Halo. 

In the next Section we test some of these interpretations 
by means of synthetic CMDs. 

6.4 Synthetic Color-Magnitude Diagrams and 
Model Luminosity Functions 

Synthetic CMDs have been generated using the TRILEGAL 
code described in Girardi et al. (2005). The detailed proce- 
dure is outlined in Carraro et al. (2002). Using typical val- 
ues derived from our observations (see Sect. 3), we have also 
simulated the photometric errors as a function of B and R 
magnitude. 



Using our results of the isochrone fitting, we generated 
a synthetic CMD for a cluster which underwent an instan- 
taneous burst of star formation 1.7 Gyr ago, and which has 
a population of 290 evolved stars. This number of evolved 
stars was derived from the luminosity function discussed 
in Sect. 4; and it is somewhat uncertain due to (1) possi- 
ble errors while removing the contamination by the field, 
and (2) Poisson statistics. If we neglect for the moment 
the uncertainties in the field contamination, and assume the 
Kroupa (2001) Inital mass Function (IMF) corrected for bi- 
naries, this number of evolved stars implies a cluster mass 
of 33.3 ± 2.0 x 10 3 M Q . 

Our initial CMD simulation includes a fraction f=30% 
of detached binaries, and assumes that their mass ratio is 
uniformly distributed between 0.7 and 1.0. This assumption 
is in agreement with the observational data for the LMC 
clusters NGC 1818 (Elson et al. 1998), NGC 1866 (Barmina 
et al. 2002), and NGC 2173 (Bertelli et al. 2003). It is worth 
recalling that the photometry of a binary with mass ratio 
smaller than 0.7 is almost indistinguishable from its primary 
alone, so that extending the interval of simulated mass ra- 
tios would not change the results. 

The results for this initial choice of parameters are 
shown in Fig. 10. The left panel shows the observational, 
decontaminated, CMD of NGC 2154 as derived in the pre- 
vious section, and the right panel shows the result of our 
best simulation; chosen among many synthetic CMDs gen- 
erated with the same input parameters but varying the ran- 
dom seed (Bertelli et al. 2003). The top panel presents the 
observational (see Sect. 4.2) ILF (Integrated Luminosity 
Function, solid line) together with the corresponding best 
fitting model ILF (dashed line). 

In the simulations, special attention was given to the 
shape of the TO region, and to the color and magnitude 
of the RGB clump. Because of the complicated structure 
of the TO, which is broadened by the presence of binaries 
and the extended star formation period, for the derivation 
of the cluster's distance and reddening we have used mainly 
the RGB clump. As can be readily seen, the mean color 
and magnitude of the clump has been well reproduced, al- 
lowing us to infer a reddening E(B — R) — 0.09 ± 0.02 
(E(B - V) = 0.057), and distance modulus (B - M B )o = 
(V - M v )o = 18.48±0.10. Both these values are well within 
the widely accepted estimates for the LMC. 

The red clump stars, however, seem to present a more 
elongated, tilted, distribution in the CMD, than in the ini- 
tial simulation. This kind of structure might be caused by 
an increased fraction of binaries, by differential reddening, 
or by an age spread inside the cluster. All of these effects 
would affect not only the CMD, but also its integrated LF, 
and in different ways. 

We have investigated these possible effects by run- 
ning additional simulations in which we varied the following 
quantities: 

• The fraction of binaries /, from (no binaries) to 1 
(each star drawn from the IMF is the primary of a system 
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with mass ratio between 0.7 and 1) at steps of 0.1. Because 
the binaries are added with equal probability to both the 
MS and red giant part of the CMD, / affects the predicted 
dwarf/giant ratio, and hence the ILF. 

• The la dispersion in reddening, between and a max- 
imum of 0.09 in E(B — R). Of course, negative values of 
reddening are not allowed, and become 0. This causes a 
maximum broadening of about 0.2 mag of the red clump 
in B — R, without affecting the overall shape of the ILF. 

• The duration of the star formation episode that gen- 
erated the cluster, from to 2.0 Gyr, with steps of about 
0.2 Gyr, centered at an age of 1.75 Gyr. This spread has a 
modest effect in the red clump luminosity, but affects the 
turn-off region of the ILF. 

For each model in the grid, we compute the reduced x 2 > de- 
fined as the mean squared difference between model and 
observations, in the magnitude interval 16 < B < 21.5. 
Fainter stars are not included both because (1) below the 
turn-off region the ILF becomes very sensitive to the IMF, 
which is not known with enough accuracy; and (2) at those 
magnitudes the completeness corrections become large and 
consequently more uncertain. 

The model with the lowest value of \ 2 (=347) has 
/ = 0.7, aE(B - R) = 0.09 mag, and a SFH duration 
of 1.2 Gyr. The corresponding best fitting synthetic CMD 
and model ILF are presented in Fig. 11. The top panel 
presents the observational (see Sect. 4.2) ILF (solid line) 
together with the model ILF (dashed line). This theoretical 
ILF therefore represents the same population, with the same 
binary fraction, same star formation history, and photomet- 
ric errors as the synthetic CMD, and with as many red giants 
as in the observations. It can be readily noticed that they 
agree reasonably well down to B — 23.5, magnitude level at 
which the incompleteness corrections are still smaller than 
50%. It should be noticed however, that the differences be- 
tween them are statistically significant, i.e. larger than the 
la error bars (67% confidence level) determined from the 
Poisson statistics. 

We would like to note that a solution of similar quality, 
with x 2 = 549, is found for / = 0.4, aE(B - R) = 0.0 mag, 
and a SFH duration of 2 Gyr. This second solution has rea- 
sonable values of / (well in agreement with estimates for 
other LMC clusters) and aE(B — R). Despite the good qual- 
ity of the ILF fit as measured by the \ 2 , a 2 Gyr duration for 
the star formation in NGC 2154 seems unacceptably long. 



7 CONCLUSIONS 

We have presented and discussed deep BR photometry of 
the intermediate-age open cluster NGC 2154 in the LMC. By 
using theoretical tools, namely isochrones, synthetic CMDs, 
and an ILF, we obtained estimates of the cluster fundamen- 
tal parameters. The distance and reddening found fall within 
the commonly accepted values for the LMC. 
Very interestingly, we found that the cluster CMD, and in 
particular the ILF, can be properly interpreted only allowing 
for a extended period of star formation (1.2 Gyr). Another, 



less extreme, example of extended period of cluster forma- 
tion, is that of the almost coeval cluster NGC 2173 in the 
LMC, with 0.3 Gyr (Bertelli et al. 2003). Whether the pres- 
ence of such extended period of star formation is physically 
possible, is to be investigated in future works. From our 
side, it is mandatory to mention that the detection of such 
a prolonged star formation period could well be caused by 
inadequacies in the stellar evolutionary tracks themselves. 
The most obvious among these inadequacies is the uncer- 
tain efficiency of core convective overshooting in MS stars 
of masses between 1 and 2 Mq (see e.g. Chiosi 2006). 

An interesting point raised by the referee is whether 
the blue stars which concur to enlarge the MS and 
thicken the TO producing the extended star forma- 
tion effect we find, might indeed be Blue Stragglers 
Stars (BSSs). These stars are ubiquitous, and they 
are routinely found in Dwarf Galaxies, Galactic star 
clusters and the general Galactic field (de Marchi 
et al. 2006). Unfortunately, a comprehensive search 
and analysis of BSSs in the LMC clusters is still un- 
available. 

While there is no consensus yet on the mechanism 
which produces these stars, their position in the 
CMD and their population are reasonably well un- 
derstood. They occupy a strip along the extension 
of the MS above the TO point, but tend to be bluer 
than the standard binary star sequence (which runs 
parallel to the MS). Besides, they occupy a region 
of the CMD where young stars of the general stel- 
lar field toward a star cluster or a dwarf galaxy are 
found. This complicates their detection and demand 
more effective membership criteria. 
In our CMD, after the cleaning procedure, most of 
blue stars are actually more compatible with field 
stars, and the shape of the TO point does not seem 
to be affected by classical BSSs, which instead would 
occupy a region which detaches from the cluster MS 
at B ~21.5, (B-R) ~0.35, following the extension of 
the cluster MS below the TO. 

The two classical explanations for these stars are 
that they are binary stars or more massive stars 
born in a separate later star formation episode. Our 
simulations include both these effects in a way that 
it is impossible to clarify whether BSSs are present 
or not. 

We have visually compared the CMD of NGC 2154 
with the CMD of the rich star cluster NGC 2173 
by Bertelli et al. (2003). The main motivation for 
this choice is that this study investigates a cluster 
-NGC 2173- which is coeval to NGC 2154 and, in- 
terestingly, the authors find that only a sizable age 
dispersion can explain the shape of the TO. 
As in our case, the cleaned CMD of this cluster 
(their Figs 4 and 5) show only a few stars on the 
blue side of the MS, and the region right above the 
TO is more easily explained in terms of binary stars 
and extended star formation. 

The fact that another coeval cluster -NGC 2173- 
does show the same features as NGC 2154 and 
can be only interpreted as having experienced long- 
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lasting star formation might instead be telling us 
that the physics of TO mass stars typical of this age 
might have some problems, as mentioned above. 
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Figure 10. The observational CMD of NGC 2154 is shown in the left panel. The right panel presents a synthetic CMD of the cluster 
obtained for our initial choice of parameters (i.e., assuming Z = 0.005, an age of 1.7 Gyr, and 30% of binaries). The top panel presents 
the observational ILF (solid line) together with best model ILF obtained with the above parameters (dashed line), la Poisson error bars 
are included. 



12 Baume et al. 



3000 



2000 



1000 



I h 



-EE ! EE— 



18 



16 



20 



22 



B 



16 



24 



18 



18 



m 20 



20 



22 



22 



24 



-10 12 3 
B-R 



24 



12 3 
B-R 



Figure 11. The observational CMD of NGC 2154 is shown in the left panel. The right panel presents our best fitting synthetic CMD 
of the cluster. The top panel presents the observational ILF (solid line) together with our finally adopted model ILF (dashed line), la 
Poisson error bars are included. See text for details. 



